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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


The  central  focus  of  the  research  performed  under  this  contract  involved 
3  1  3+ 

the  optical  transition  of  the  Pr  ion  in  LaF^.  Our  plan  was  to 

perform  a  photon  echo  modulation  experiment  similar  to  the  one  used  to  study 

3  3  3+1 

the  optical  -  Pq  transition  of  the  Pr  ion  in  LaF^.  (In  that 

experiment  modulation  data  provided  ground  and  excited  state  spectroscopic 

information,  nuclear  hyperfine  linewidths,  the  optical  homogeneous  linewidth 

and  the  relative  orientation  of  the  principle  axes  of  the  ground  and  excited 

state  hyperfine  hamlltonlans. 

3+ 

The  Pr  ion  in  crystals  of  LaF^  substitutes  for  the  La  atom  and 

exhibits  an  absorption  and  emission  spectrum  characteristic  of  an  ion  whose 

2 

energy  levels  are  nondegenerate.  Although  the  crystal  structure  of  LaF^  is 

somewhat  controversial,  the  most  likely  model  appears  to  be  that  proposed  by 

3  4  3 

Mansmann  and  Zalkln  et  al.  They  suggest  P3cl  -  which  has  a  trigonal 

structure  with  six  molecules  per  unit  cell  and  which  exhibits  C2  point 

symmetry  about  each  of  the  three  inequivalent  La  sites.  The  twofold  axis  of 

Cg  is  normal  to  the  crystal  c  axis.  Since  the  Pr  ion  has  two  unpaired  4f 

electrons,  the  low  symmetry  at  the  La  site  removes  all  the  degeneracy  of  the 

energy  eigenstates  and  requires  that  for  each  eigenstate  <J^>«  5  where  J  is 

3+ 

the  electronic  angular  momentum  operator.  To  first  order  then,  the  Pr  ions 
have  no  electronic  magnetism  and  one  might  expect  that  they  would  be  well 
Isolated.  It  follows  that  in  any  coherent  optical  transient  experiment  they 
should  behave  in  a  simple  manner.  It  is  not  surprising  then  that  Takeuchl 

C  £ 

and  Szabo  found  that,  in  contrast  to  the  case  of  ruby,  it  is  possible  to 

u 

observe  photon  echoes  in  LaF^:Pr  without  the  application  of  an  external 

3  3  0 

magnetic  field.  They  worked  in  the  -  Pq  transtiion  (4777A)  at  liquid 


htlltm  temperatures.  The  echoes  appeared  to  decay  exponentially  with  pulse 
separation  t. 

The  first  observation  of  a  nonexponential  decay  for  the  echo  intensity 
3+  7 

versus  t  in  LaF^:Pr  was  reported  by  Chen  et  al.  who  found  a  sharp 

rephasing  at  t*  240  nsec,  just  beyond  the  region  studied  by  Takeuchi.^  They 

13  0 

also  reported  the  observation  of  echoes  In  the  transition  (5925  A) 

3 

which  were  modulated  at  the  inverse  of  8.47  MHz,  the  ground-state 

O 

splitting  frequency.  More  recent  experiments  by  Chen  et  al.  and  Morsink  and 

9 

Vlersma  using  independently  triggered  nitrogen-pumped  dye  lasers  have 

confined  the  observed  rephasing  at  240  nsec  and  have  shown  that  the 

echo-decay  envelope  is  more  complex  than  had  previously  been  expected.  This 

complexity  is  fortunate  however  as  it  allows  analysis  to  discover  the  most 

u 

complete  description  of  the  Pr  ion  with  its  environment. 

Our  first  work  on  the  ^H^  -  *D2  transition  encountered  difficulty 
because  of  the  relatively  weak  character  of  the  transition  dipole  moment. 
Experimental  difficulty  aside  we  were  able  to  obtain  echo  modulation  data  but 
we  were  unable  to  analyze  it  successfully.  Two  problems  were  encountered. 
Plrst  the  best  theoretical  fit  we  could  obtain  to  the  modulation  data  was  not 
very  satisfying  and  this  fit  required  that  we  discount  Erickson's^’ ** 
measurements  of  the  *D2  state  splitting.  On  re-analyzing  Erickson's  work  we 
felt  (as  did  Erickson  in  a  private  conversation)  that  for  the  *D2  state  the 
results  were  not  unambiguous  and  a  more  direct  measurement  of  the  excited 
state  splitting  was  in  order.  In  Erickson's  experiment  he  ’‘sees”  many 
resonances  and  these  must  be  carefully  handled  to  achieve  the  desired  result. 
There  was  also  the  surprising  circumstance  that  his  inferred  excited  state 
splittings  of  3.7  MHz  and  4.7  MHz  summed  too  close  to  the  8.48  MHz  splitting 
of  two  of  the  ground  state  levels.  This  coincidence  was  disturbing  to  both 
Erickson  and  us.  We  therefore  proceeded  to  make  a  direct  measurement  of  the 


terminal  level  apllttinga  by  the  PENDOR  technique.  This  method  was  developed 

in  our  laboratory  and  was  successfully  used  to  obtain  the  hyperflne 

3+ 

interactions  between  the  Cr  ion  and  the  neighboring  A1  nuclei  in  ruby 
3+ 

(Cr  :A1203).  The  method  has  the  advantage  that  its  results  are  unambiguous. 
It  works  by  exploiting  the  dependence  of  the  three-pulse  stimulated  photon 
echo  amplitude  on  the  application  of  an  rf  field  which  is  resonant  with  a 
hyperflne  transition.  We  performed  such  an  experiment  and  found  Erickson's 
work  to  be  wholly  correct.  This  work  was  published  in  Physical  Review  and 
appears  herein  as  Appendix  A. 

Our  difficulty  in  analysing  our  photon  echo  modulation  data  was  traced 

to  the  existence  of  a  long-lived  stimulated  echo  which  constructively 

interfered  with  our  two-pulse  photon  echo.  In  our  experiment,  the  ordinary 

two  pulse  echo,  needed  to  gather  the  modulation  data,  is  formed  by  a  pair  of 

pulses,  separated  in  time  by  0.1  microseconds  to  9  microseconds.  The 

effect  of  the  two  pulses  is  to  produce  the  ordinary  photon  echo  delayed  from 

the  second  pulse  by  the  pulse  separation.  The  pulse  pair  also  induces  order 

in  the  population  distribution  of  the  ground  state  hyperflne  levels  (optical 

pumping)  and  due  to  the  very  long  spin  lattice  relaxation  time  of  the 
12 

system,  the  order  will  persist  for  as  long  as  several  minutes.  A  laser 
pulse  applied  to  such  an  ordered  system  will  generate  a  stimulated  echo, 
delayed  in  time  from  the  pulse  by  the  same  amount  as  the  pulse  separation 
between  the  pair  of  pulses  which  induced  the  order  in  the  first  place.  In 
our  usual  data  taking  method,  pairs  of  pulses  arrive  at  the  sample  at  the 
repetition  frequency  of  the  laser  (10  Hz.),  forming  a  complicated  sequence  of 
echoes  and  stimulated  echoes  simultaneously.  To  eliminate  this  problem  we 
applied  radio-frequency  energy  at  one  of  the  ground  state  hyperflne 
resonances  (16.7  MHz)  during  the  interval  between  laser  shots,  effectively 
eliminating  the  population  distribution  order  and  the  long-lived  stimulated 


Eliminating  the  stimulated  echo  removed  the  major  obstacle  in  the  way  of 

obtaining  reliable  modulation  data.  The  other  major  problem  area  was  signal 

to  noise.  Echo  modulation  experiments  require  signal  handling  with  a  wide 

4 

dynamic  range.  He  were  ultimately  able  to  see  echoes  over  a  range  of  10  but 

that  was  only  possible  through  hard  work.  Narrowing  the  spectral  width  of 

our  laser  to  4  GHz  helped  as  well  as  going  to  samples  whose  resonances  were 

just  on  the  verge  of  being  optically  thick.  We  Improved  our  computer 

assisted  data  taking  process.  In  particular,  two  new  computer  operated 

peripherals  were  Installed  in  the  system.  The  first  Is  a  computer 

controllable  electrical  attenuator  which  Is  used  to  maintain  the  modulated 

signal  within  the  linear  range  of  the  data  acquisition  electronics.  The 

computer  Is  thus  able  to  adjust  the  signal  level  automatically.  The  second 

device  is  a  pair  of  digitally  programmable  pulse  generators,  used  to  set  the 

laser  pulse  separation  and  Pockel's  cell  gate.  These  devices  are  of  a 

“random  access**  nature  allowing  the  pulse  separation  to  be  varied  over 

several  orders  of  magnitude  between  a  single  laser  shot.  This  allows  us  to 

periodically  measure  the  echo  magnitude  at  a  single  fixed  reference  time 

without  suffering  a  loss  In  the  data  taking  rate.  The  reference  measurement 

allows  us  to  compensate  for  echo  amplitude  drift  during  the  course  of  a 

single  data  taking  run,  as  well  as  providing  a  normalisation  factor  which 

facilitates  averaging  data  sets  from  different  runs. 

With  the  Improvements  In  the  technique  outlined  above  we  have  been  able 

3  1 

to  obtain  excellent  echo  modulation  data  with  which  to  analyze  the 

3+ 

transition  of  P r  In  LaF^.  The  analysis  involves  both  Fourier  analysis  to 
yield  spectroscopic  information  directly  and  theoretical  echo  modulation 
reconstruction  to  provide  the  detailed  form  of  the  interaction  hamiltonian. 
This  work  will  be  published  In  Physical  Review  and  a  proof  of  the  article 


appears  la  Appendix  B. 

The  resonance  widths  of  Pr^sLaF^  reflect  the  total  (lnhoeogeneous  and 

homogeneous)  nuclear  line widths  of  the  various  hyperfine  transitions. 

1  13 

Previous  measurements  '  have  shown  the  ground  state  widths  to  be  abo.it  200 

kHs.  This  width  is  wostly  lnhoeogeneous  and  arises  froe  the  interaction  of 

141 

the  enhanced  ground  state  Pr  nuclear  noaent  of  11.3  kHz/G  with  the 

19 

surrounding  distribution  of  static  P  nuclear  eooents.  Mac far lane  and 
ShcJ-hy1^  have  aeasured  the  g-tensor  for  the  lD2  state  and  found  the 
enhanceaent  to  be  roughly  one-fifth  as  much  as  in  the  ground  state.  This 
would  lead  one  to  expect  a  auch  narrower  llnewidth  for  the  nuclear 
transitions  in  the  Dj  state  and  our  measurement,  while  very  nearly 
resolut ion-1 lal ted  by  the  overall  echo  decay  envelope,  establishes  a  new 
upper  Halt  of  30  kHz  [FWHM]  for  the  hyperf ine  line  at  8.51  Mhz,  consistent 
with  this  expectation.  The  llnewidths  of  the  hyperfine  at  3.72  and  4.79  MHz 
are  somewhat  larger,  being  70  kHz,  and  60  kHz  respectively.  The  results  for 
the  lD2  excited  state  paraaeters  are  suaaarized  in  Table  1. 

The  effective  hyperfine  Haslltonian  is  given  by 

Hg(e)  "  Pg(e)^Xz(Z)+  ~S  (Ix(X)  ”  Xy(Y) ^ 

where  Y<«>  •nd  sre  Interaction  constants.  The  upper  and  lower  case 

subscripts  eaphaslze  that  the  ground  (g)  and  excited  (e)  state  principal  axes 
need  not  be  the  saae,  although  Neumann' s  principle  requires  that  they  have  at 
least  one  coamon  axis  parallel  to  the  local  C2  syaaetry  axis.  This 
constraint  significantly  reduces  the  number  of  possible  orientations  and  we 
••y  represent  these  possibilities  on  a  two-dimensional  map.  The  coordinates 
of  such  a  map  give  the  particular  choice  of  common  axis,  and  the  angle  9 
through  which  one  set  of  axes  is  rotated  with  respect  to  the  other  set  about 


the  cowon  axis*  Since  the  Haailtonlan  in  Eq.  1  is  Insensitive  to  the  sign 
of  the  coordinates  an  orientation  of,  say,  x  ||  X  is  equivalent  to  x | | -X.  We 
have  calculated  the  echo  modulation  in  increments  of  five  degrees  in  the 
rotation  angle  8.  In  the  sup  presented  in  Pig*  1,  for  each  orientation 
coordinate  a  box  is  drawn  whose  width  is  proportional  to  the  minimum  mean 
square  deviation  between  our  dita  and  the  calculated  modulation*  Only  the 
first  microsecond  of  data  were  used  in  the  fitting  procedure  for  Figure  1. 
Beyond  one  ysec.  the  data  are  dominated  by  the  excited  state  hyperfine 
splittings,  due  to  the  relative  llnewidths  of  ground  and  excited  state 
levels,  and  Including  these  data  in  the  fit  considerably  reduces  the 
selectivity  of  the  map. 

Figure  2  shows  in  expanded  scale  the  first  three  psec  of  the  data,  and 
Figure  3  shows  the  theoretical  curve  that  was  calculated  using  the 
orientation  coordinate  giving  the  best  fit  to  the  data  (as  indicated  by  an 
arrow  in  the  lower  left  of  Figure  i). 

The  map  in  Figure  1  shows  a  high  degree  oi  symmetry  and  concommitant 
lack  of  selectivity  for  determining  the  orientstion  parameter.  To  explore 
the  extent  to  which  this  symmetry  is  intrinsic,  we  have  substituted  for  the 
data  used  to  calculate  the  map  of  Figure  1  the  values  of  the  theoretical 
modulation  corresponding  to  Figure  3,  producing  the  map  shown  in  Figure  4. 

The  scale  in  Figure  4  was  chosen  to  make  the  box  corresponding  to  the 
coordinate  with  the  worst  fit  in  Figure  4  approximately  the  same  size  as  the 
box  with  the  worst  fit  in  Figure  2.  It  is  apparent  that  the  high  degree  of 
symmetry  persists.  To  estimate  the  extent  to  which  our  signal-to-noise  would 
have  to  improve  to  allow  us  to  select  one  orientation  with  certainty,  we 
replotted  Figure  4  with  the  box  width  scale  expanded  by  a  factor  of  ten,  as 
shown  in  Figure  5.  It  appears  that  an  improvement  in  signal-to-noise  by  a 
factor  of  ten  or  twenty  would  have  made  possible  a  more  definitive  statement 


•bout  Che  orientation.  This  improvement  will  probably  come  from  the  use  of 
frequency  and  amplitude  stabilised  dye  lasers.  Frequency  jitter  in 
particular  causes  considerable  fluctuation  in  echo  intensity,  presently 
averaged  by  the  data^taklng  system,  and  such  averaging  tends  to  wash  out  some 
of  the  fine  details  of  the  modulation  pattern.  Despite  the  high  symmetry,  we 
can  draw  a  number  of  conclusions  about  the  relative  orientation.  All 
configurations  with  all  axes  parallel  are  ruled  out,  as  are  those  with  the 

excited  state  Z-axls  parallel  to  C2*  These  conclusions  are  similar  to  the 

1  3  3 

measurement  of  Chen  et^  al. ,  on  the  relative  orientation  of  the  and  Pq 

1  13 

states.  Since  the  g-tensor  is  nearly  isotropic  in  the  state  and  absent 
3 

completely  in  the  Pq  state  the  relative  orientation  parameter  must  be 

primarily  determined  by  the  relative  orientation  of  the  pure  quadrupole 

tensor  with  respect  to  the  ground  state  hyperfine  tensor. 

3  13+ 

In  addition  to  the  study  of  the  -  D2  Pr  transition  we  also 
performed  a  theoretical  analysis  of  propagation  narrowing  in  the  transmission 
of  a  light  pulse  through  a  spectral  hole.  This  work  appeared  in  Physical 
Review  and  is  attached  as  Appendix  C.  Because  of  the  linearity  of  the 
aystem,  it  is  possible  to  treat  the  pulse  as  if  it  were  Interacting  with  two 
entirely  separate  groups  of  atoms.  The  first  group  comprises  the  usual  dipole 
oscillators  in  the  full  Doppler  line,  and  the  second  group  is  a  set  of 
"negative"  oscillators  in  a  narrow  region  which  serve  to  "eat  out"  the  hole. 
The  negative  oscillators  emit  rather  than  absorb  light.  It  is  the  emission 
of  the  "negative"  oscillators  that  cause  narrowing  of  the  transmitted  pulse 
spectrum  as  obtained  by  our  theoretical  analysis. 
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Photon-echo  nuclear  double  resonance  in  LaFj:Pr3+ 
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Using  the  photon-echo-nuclear-double-resonance  technique,  we  have  measured  the  nuclear 
hyperfine  splittings  in  the  excited  'Dj  slate  of  LaF3:Pr3+.  Our  measurements  are  in  agreement 
with  the  results  obtained  by  Erickson  from  his  saturated-and-enhanced-absorption  study  of  the 
iHA-'D2  transition  in  LaF3:Pf34\ 


Photon-echo  modulation  in  LaF3:Pr3+  has  been  a 
valuable  source  for  obtaining  the  various  spectroscop¬ 
ic  and  relaxation  parameters  associated  with  the 
3//4 -J/>o  transition.1  This  technique  when  applied  to 
the  >Ht-'D2  transition,  however,  has  only  had  modest 
success.2  In  particular,  echo  modulation  data  taken 
in  the  standard  manner  are  not  altogether  consistent3 
with  the  measured  excited-state  splittings  as  obtained 
by  Erickson.4-5  It  was  therefore  felt  that  an  indepen¬ 
dent  measurement  of  the  'D2  excited-state  splittings 
was  necessary  to  build  confidence  in  their  accuracy. 

To  this  end  we  have  performed  a  photon-echo- 
nuclear-double-resonance  (PENDOR)  experiment6 
on  the  }Ht-'D2  transition  of  this  system. 

The  PENDOR  technique  was  first  applied  to  the 
measurement  of  energy  splittings  of  I3A1  nuclei  in  a 
single  crystal  of  ruby,  and  was  used  to  obtain  the  in¬ 
teraction  parameters  of  the  Cr3+  ion  with  its  sur¬ 
rounding  nuclear^neighbors  in  both  the  ground  2Aa 
and  excited  2E(E)  states  of  ruby.7  Until  our  use  of 
this  technique  to  measure  nuclear  hyperfine  splittings 
in  LaF3:Pr3+,  PENDOR  had  not  been  applied  to  the 
study  of  any  material  other  than  ruby. 

A  PENDOR  experiment  may  be  summarized  as 
follows:  Three  laser  pulses  of  resonant  radiation  are 
made  incident  on  a  sample  to  produce  a  stimulated 
photon  echo,  with  a  pulsed  rf  field  applied  to  the 
sample  during  the  interval  between  the  second  and 
third  pulses.  When  the  rf  field  is  resonant  with  a  nu¬ 
clear  transition  in  either  of  the  two  terminal  levels  of 
the  echo  transition,  there  will  be  a  degradation  of  the 
echo  intensity.  This  decrease  in  echo  size  is  the  dou¬ 
ble  resonance  signal.  As  is  generally  true  in  optical-rf 
double-resonance  experiments,  one  obtains  a  sensi¬ 
tivity  greatly  increased  over  purely  NMR  or  optical 
methods  for  several  reasons.  With  PENDOR.  the 
resolution  is  limited  not  by  the  linewidth  of  the  laser 
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used,  as  was  the  case  in  Erickson’s  enhanced-and- 
saturated-absorption  method,4  but  only  by  the  length 
of  time  over  which  the  rf  field  is  applied  to  the  sam¬ 
ple.  Secondly,  there  are  no  ambiguities  in  the  inter¬ 
pretation  of  the  observed  line  shape  For  example,  in 
order  to  obtain  the  true  inhomogeneous  line  shape  of 
the  nuclear  transitions  in  the  'D2  state  in  LaF3:Pr3+, 
Erickson  had  to  deconvolve  the  measured  line  shape 
from  the  linewidth  of  the  laser  that  he  used.  The  in¬ 
ferred  linewidths  are  quite  sensitive  to  whether  the 
measured  line  shape  is  assumed  to  be  Gaussian  or 
Lorentzian.  A  third  and  important  advantage  of 
PENDOR  is  that  it  can  be  used  to  study  nuclear  hy¬ 
perfine  structure  in  excited  states  of  ions  in  solids 
that  are  not  accessible  to  cw  dye  lasers  because  of 
power  or  spectral  range  considerations. 

Stimulated  photon  echoes  were  produced  on  the 
transition  (592.5  nm)  connecting  the  lowest  levels  of 
the  crystal-field-split  components  of  the  3 //4  and  'D2 
states.  These  levels  are  electronic  singlets  that  are 
further  split  into  three  doubly  degenerate  nuclear  lev¬ 
els  by  the  combined  influence  of  the  Pr3+  nuclear 
quadrupolar  and  second-order  magnetic  hyperfine  di- 
po'  r  interactions.  The  previously  measured1-4  nu¬ 
clear  splittings  for  both  the  3 //4  and  ' D2  states  are 
small  compared  to  the  10-GHz  bandwidth  of  the  laser 
excitation  pulses.  The  timing  sequence  for  a  PEN¬ 
DOR  experiment  is  shown  in  Fig.  1. 

Our  experimental  setup  utilized  two  independently 
triggered  nitrogen-laser-pumped  dye  lasers,  which 
produce  multiplo-kW,  5-nsec-long  pulses  in  a 
bandwidth  of  approximately  10  GHz.  The  output  of 
one  of  the  lasers  was  divided  into  two  beams  of 
roughly  equal  intensity  by  beam  splitter  BS  (see  Fig. 
2).  One  beam  was  steered  into  an  optical  delay  line 
(consisting  of  a  White  cell  with  spherical  mirrors  hav¬ 
ing  radius  of  curvature  equal  to  2.9  m).  while  the 
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FIG.  1.  Laser  and  radio-frequency  excitation  pulses  for 
PENDOR  experiment. 

other  was  combined  and  made  collinear  with  the  out¬ 
put  of  the  second  laser  by  beam  combiner  BC1.  The 
delayed  beam  from  the  White  cell  was  then  combined 
with  the  other  two  beams  at  BC2.  The  temporal 
separation  between  laser  pulses  1  and  2  was  adjusted 
to  be  130  nsec,  and  the  separation  between  pulses  2 
and  3  was  adjusted  to  be  either  20  or  40  fistc.  The 
three  collinear  pulses  were  focused  to  a  100- pi  m- 
diameter  spot  in  the  sample  by  a  20-cm-focal-length 
lens.  The  sample,  whose  PrJ+  concentration  was  1 .0 
at.  wt.%,  had  a  thickness  of  2.5  mm,  and  was  mount¬ 
ed  in  a  cryostat  that  was  maintained  at  a  temperature 
of  2.5  K.  The  c  axis  of  the  crystal  was  perpendicular 
to  its  optically  buffed  surface,  and  parallel  to  the 
direction  of  propagation  of  the  laser  pulses.  The  rf 
pulse  was  applied  to  the  sample  by  means  of  a  1-/aH, 
4-mm-diameter  coil  consisting  of  10  turns  of  30- 
gauge  copper  wire  wound  about  (he  sample  in  such  a 
way  that  the  oscillating  rf  field  was  parallel  to  the 
crystal  c  axis.  The  coil  was  energized  by  a  broadband 
3-W  amplifier  (ENI  model  300  L)  driven  by  a 
Hewlett-Packard  frequency  synthesizer  whose  output 
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FIG.  2.  Experimental  arrangement  for  PENDOR  mea¬ 
surements.  See  text  for  detailed  explanation. 


was  gated  on  between  the  second  and  third  laser 
pulses.  The  rf  gating  circuit  was  triggered  by  a 
multipulse-generator  unit  that  was  interfaced  to  a 
PDP  8/e  minicomputer,  making  it  possible  to  gate 
the  rf  field  on  and  off  under  program  control.  The 
output  of  the  frequency  synthesizer  was  also  con¬ 
trolled  by  the  PDP  8/e.  To  enhance  the  rf  field  at 
the  sample,  an  impedance  matching  circuit  was  in¬ 
serted  between  the  rf  amplifier  and  the  coil  at  the 
sample.  After  passing  through  the  sample,  the  laser 
beams  were  recollimated  by  another  20-cm-focal- 
length  lens  and  then  directed  into  an  RCA-C3I034 
photomultiplier  tube  (PMT).  The  three  excitation 
pulses  were  prevented  from  saturating  the  PMT  by 
three  stages  of  crossed  polarizers  and  Pockels  cells 
that  were  pulsed  open  about  30  nsec  before  the  ar¬ 
rival  time  of  the  stimulated  echo.  The  output  of  the 
PMT  was  amplified  and  fed  into  a  gated  stretcher  and 
then  directed  to  an  analog-to-digital  converter.  To 
eliminate  the  effects  of  long-term  drift  in  the  echo 
signal  due  to  changes  in  laser  intensity  or  optical 
alignment,  the  rf  field  was  alternately  gated  on  and 
off,  with  the  computer  separately  averaging  the 
stimulated  photon  echo  intensity  for  the  two  cases. 
The  average  was  taken  over  200  shots.  The  value  of 
(/.i  -/nrf)//nrf.  where  / rf ( f nrf )  stands  for  the  echo  in¬ 
tensity  with  (without)  the  rf  field  gated  on,  is  then 
printed  out  by  the  computer  as  a  function  of  rf  fre¬ 
quency.  The  Q  of  our  matching  circuit  was  ~  100 
and  we  had  to  adjust  the  matching  circuit  capacitors 
whenever  the  frequency  was  changed. 

Figure  3  shows  the  percentage  decrease  in  the 
stimulated  photon  echo  intensity  as  a  function  of  the 
applied  rf  frequency.  The  dips  at  3.7  and  4.65  MHz 
correspond  to  the  energy  splittings  in  the  excited  '£>2 
state,  and  agree  well  with  Erickson’s  reported  values 
of  3.7  and  4.7  MHz.4  Our  measured  linewidths  of 
350  kHz  (full  width  at  half  maximum  (FWHM)]  for 
the  two  resonances  are,  however,  larger  than  the  200 
±  50-kHz  linewidth  that  he  reported.8  This  discrep- 
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FIG.  3.  PENDOR  resonance  curve,  showing  ’/>2  exciled 
stale  splittings  in  LaF3:PrJ+ 
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ancy  can  perhaps  be  explained  by  the  suppositon  that 
the  200-kHz  linewidth  he  inferred  is  a  consequence 
of  his  assuming  a  Lorentzian  line  shape  for  decon¬ 
volving  the  observed  line  shape  from  the  laser  line 
shape.  For  a  Gaussian  line  shape,  the  inferred 
linewidth  can  be  a  factor  of  2  larger.9  In  Fig.  3  the 
solid  and  dotted  lines  are  Gaussian  and  Lorentzian 
fits,  respectively,  to  the  data  points  with  a  (FWHM) 
linewidth  of  350  kHz.  The  quality  of  our  data 
prevents  our  preferring  one  over  the  other,  but  since 
a  large  part  of  the  inhomogeneity  of  the  nuclear  hy- 
perfine  levels  is  due  to  the  perturbing  effects  of 
eleven  19F  neighbors,10  one  would  expect  the 
linewidth  to  be  Gaussian  rather  than  Lorentzian. 

Previous  work  which  Fourier  analyzed  photon- 
echo-modulation  data  for  the  }HA-iP0  transition1  has 
shown  that  the  inhomogeneous  linewidth  of  the  nu¬ 
clear  transitions  in  the  1P<>  state  is  approximately  20 
kHz,  while  the  inhomogenous  linewidth  of  the  3HA 
ground  state  is  approximately  200  kHz.  For  the  }Po 
state,  7  —  0,  and  the  linewidth  is  due  primarily  to  the 
static  dipolar  and  quadrupolar  interactions  of  the  l4,Pr 
nucleus  with  the  neighboring  ,9F  nuclei  and  the  elec¬ 
tric  field  gradient  at  the  l41Pr  site,  respectively.  For 
the  3H 4  state,  7—4  and  the  close  spacing  of  the 
crystal-field-split  levels  of  the  3HA  multiplet  gives  rise 
to  an  enhanced  nuclear  moment  that  interacts  strong¬ 
ly  with  the  surrounding  nuclear  neighbors  and  leads 
to  an  increased  linewidth.  In  like  manner  the  lD2 
multiplet  levels  are  broadened  beyond  the  20-kHz 
width  of  the  3P0  state.  According  to  a  second  mo¬ 


ment  calculation,  the  contribution  to  the  inhomo¬ 
geneous  linewidth  from  the  unenhanced  Pr-F  nuclear 
dipolar  interaction  is  13  kHz.11  When  we  take  into 
account  the  experimentally  determined  enhanced  nu¬ 
clear  gyromagnetic  ratio  of  y/2 n  -  10  kHz/G  for  the 
'D2  state1112  we  obtain  an  inhomogeneous  linewidth 
of  approximately  120  kHz.  The  contribution  from 
other  sources  such  as  lifetime  broadening  ( —  600 
Hz)  and  phonon  processes  ( —  800  Hz  )  (Ref.  4)  are 
negligible.  Since  we  typically  used  rf  pulses  20-jzsec 
long,  the  contribution  to  the  linewidth  due  to  the  fi¬ 
nite  duration  of  the  rf  field  is  8  kHz.  The  remaining 
220  kHz  of  the  measured  350-kHz  linewidth  that 
must  be  accounted  for  is  possibly  due  to  rf  power 
broadening. 

In  conclusion,  using  a  photon-echo  nuclear- 
double-resonance  technique,  we  have  shown  that  the 
nuclear  splittings  in  the  lD2  state,  as  inferred  by  Er¬ 
ickson  from  his  saturated-and-enhanced-absorption 
study  of  the  3HA-'D2  transition  in  LaFt:PrJ+,  are 
correct. 
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We  have  performed  a  photon  echo  modulation  measurement  of  the  xD2-iH4  levels  of 
Pr,+:LaFj.  This  time-domain  method  yields  a  precise  determination  of  the  [D 2  excited- 
state  hyperfine  splittings  and  linewidths.  New  information  on  the  relative  orientation  of 
the  principal  axes  of  the  excited-  and  ground-state  hyperfine  Hamiltonians  is  presented. 


In  this  paper  we  present  the  results  of  a  photon 
echo  modulation  measurement  of  the  hyperfine 
structure  and  line-broadening  character  of  the 
lowest  crystal-field-split  levels  of  the  iZ>j-3//4  tran¬ 
sition  in  Pr3+:LaF3.  These  electronic  singlets  are 
each  split  into  three  pairs  of  doublets  by  the 
second-order  hyperfine  and  electric  quadrupole  in¬ 
teractions.' 1 2  The  splittings  are  of  the  order  of  10 
MHz.  Both  the  nuclear  and  the  optical  transitions 
of  this  system  are  broadened  by  the  magnetic  in¬ 
teraction  of  the  Pr  and  F  nuclei.3  While  the  nuclear 
magnetic  interactions  are  enhanced  by  second-order 
electronic  coupling,4  these  interactions  are  still 
weak,  and  the  optical  homogeneous  linewidth  is  one 
of  the  narrowest  ever  seen  in  a  solid.5  The  spectra 
described  in  this  paper  have  been  measured  previ¬ 
ously  by  enhanced  and  saturated  absorption,2  pho¬ 
ton  echo  nuclear  double  resonance,6  and  nonmodu- 


lated  photon  echo.5  Our  result,  however,  represents 
the  first  measurement  of  the  splittings  and 
linewidths  of  this  transition  carried  out  entirely  in 
the  time  domain  and  provides  new  information  on 
the  relative  orientation  of  the  principal  axes  systems 
of  the  ground  3//4  and  excited  ]D2  hyperfine  Ham¬ 
iltonians. 

Hyperfine  spectroscopy  using  photon  echo  modu¬ 
lation  has  been  described  completely  by  Chen  et  al.1 
who  studied  the  1P0-iH4  transition  in  Pr3+:LaF3. 
We  summarize  here  the  basic  method  and  some  im¬ 
provements  we  have  made.  Two  laser  pulses, 
resonant  with  the  'Z>2-3/f4  transition  and  separated 
in  time  by  r,  are  incident  upon  a  10-mm-long  sam¬ 
ple,  cooled  to  2.5  K.  This  preparation  gives  rise  to 
a  rephased  dipole  moment  at  time  2r  after  the  first 
pulse,  given  by 


J 


-  ,  2/  +  I 

/*(2r)=P(0)(2/  -f  1 )~  j  (^W^~V^V^'lUexp[-(ft£+n&,)r] 
o,Rr,««  t 


Xcos[(&>Qy — )r ]exp(  —  2r/7’j )  , 
I - 


where  the  notation  is  that  of  Ref.  7.  The  essential 
features  of  this  formula  are  its  dependence  on  the 
W  matrices,  which  carry  information  on  the  relative 
orientation  of  the  ground-  and  excited-state  hyper¬ 
fine  Hamiltonians,  the  nuclear  hyperfine  transition 
linewidths  ft,  and  the  nuclear  hyperfine  transition 
frequencies  to.  The  measured  echo  intensity  is  the 
square  of  this  amplitude  and  will  contain  beats 
from  both  the  fundamental  hyperfine  frequencies 
and  their  harmonics. 

The  experimental  method  utilizes  two  electroni¬ 
cally  triggered  nitrogen  laser-pumped  dye  lasers. 
The  1-GHz  wide,  kilowatt  pulses  are  made  collinear 
and  focused  to  a  lOO-jxm  spot  on  the  0.1  wt.  %  sa ni¬ 


di 


pie  of  Pr^iLaFj.  The  echo  signal  is  detected 
through  a  Pockets  cell  shutter  by  an  RCA  C31034 
photomultiplier  tube.  A  PDP8/e  computer  is  used 
to  control  and  sweep  the  laser  pulse  separation 
through  a  desired  range,  and  to  integrate  and  digi¬ 
tize  the  echo  intensity  for  each  laser  shot.  To  elim¬ 
inate  the  effects  of  a  remnant  (10— 20)-nsec  laser 
pulse  separation  jitter,  a  computer-read  ORTEC 
457  time-to-pulse-height  converter  is  used  to  put 
each  echo  shot  into  the  appropriate  5-nsec-wide  bin 
in  the  computer  memory.  To  compensate  for  laser 
amplitude  drift  and  to  facilitate  averaging  data 
from  many  sweeps,  the  computer  was  made  to 
periodically  record  the  echo  intensity  at  a  fixed 
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PULSE  SEPARATION  (u»«c) 

FIG.  1.  Photon  echo  intensity  vs  laser  pulse  separation  for  the  'D^H*  transition  of  0.1  wt.  %  PrJ+:LaFJ. 


reference  time  (r =600  nsec).  Figure  1  shows  the 
average  echo  intensity  versus  pulse  separation. 
Each  point  is  the  average  of  at  least  100  shots  and 
the  data  shown  represent  the  average  of  about  30 
separate  runs.  A  given  run  covers  about  1  fisec  of 
data.  To  eliminate  the  effect  of  optical  pumping  in 
the  ground  state1  (which  for  this  method  produces 
an  artifact  known  as  long-lived  stimulated  echo8), 
the  16.7-MHz  ground-state  hyperfine  resonance  is 
saturated  with  rf  energy  between  laser  shots. 

The  overall  decay  of  the  echo  in  Fig.  1  reflects  a 
homogeneous  optical  linewidth  of  72  kHz  (corre¬ 
sponding  to  Tj=4.4  fisec  for  a  Lorentzian  line 
shape),  consistent  with  several  other  recent  measure¬ 
ments.5'’  More  information  may  be  extracted  from 
this  data  by  examining  its  Fourier  transform,  shown 
in  Fig.  2.  Before  Fourier  transformation,  the  data 
of  Fig.  1  were  multiplied  by  a  factor  exp(-f  4r/r2) 
so  that  the  resonance  linewidths  of  Fig.  2  do  not  in¬ 
clude  the  overall  optical  decay  linewidth.  The  three 
strong  resonances  at  3.72,  4.79,  and  8.51  MHz  are 
the  hyperfine  frequencies  of  the  'D2  excited  states. 
Weaker  resonances  (note  that  beyond  10  MHz  the 
vertical  scale  is  100  times  more  sensitive)  at  12.24, 
13.26,  and  17.04  MHz  are  from  mixing  of  the  first 
three  strong  resonances.  There  are  also  ground- 
state  resonances  at  8.48,  16.7,  and  25.14  MHz  but 


the  8.48-MHz  line  is  masked  by  the  stronger  8.51- 
MHz  excited-state  line,  and  the  16.7-MHz  line  ap¬ 
pears  to  be  split.  We  believe  that  this  splitting  is  an 
unanticipated  effect  of  the  16.7-MHz  rf  energy  used 
to  eliminate  long-lived  stimulated  echo,  and  prob¬ 
ably  arises  because  for  pulse  separations  greater 
than  3  fisec  the  rf  energy  was  inadvertently  allowed 
to  overlap  the  echo  pulse  sequence. 

The  resonance  widths  reflect  the  total  (inhomo¬ 
geneous  and  homogeneous)  nuclear  linewidths  of 
the  various  hyperfine  transitions.  Previous  mea¬ 
surements1,7  have  shown  the  ground-state  widths  to 
be  about  200  kHz.  This  width  is  mostly  inhomo¬ 
geneous  and  arises  from  the  interaction  of  the 
enhanced  ground-state  mPr  nuclear  moment  of 
11.5  kHz/G  (Ref.  1)  and  the  surrounding  distribu¬ 
tion  of  static  l9F  nuclear  moments.  Macfariane  and 
Shelby10  have  measured  the  g  tensor  for  the  lD2 
state  and  found  the  enhancement  to  be  roughly  5 
times  less  than  in  the  ground  state.  This  would  lead 
one  to  expect  a  much  narrower  linewidth  for  the 
nuclear  transitions  in  the  ]D2  state  and  our  mea¬ 
surement,  while  very  nearly  resolution  hmited  by 
the  overall  pcho-decay  envelope,  places  a  new  upper 
limit  of  30  kHz  (full  width  at  half  maximum)  for 
the  8. 51 -MHz  linewidth,  consistent  with  this  expec¬ 
tation.  The  3.72-  and  4.79-MHz  linewidths  are 
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FIO.  2.  Squared  magnitude  of  the  Fourier  transform  of  the  data  in  Fig.  1.  The  data  were  first  multiplied  by 
exp(r/1.2  fisec).  The  vertical  scale  is  100  times  more  sensitive  for  frequencies  beyond  10  MHz. 
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FIG.  3.  Map  of  minimum  mean-square  deviation  be¬ 
tween  echo  data  and  theoretical  curve  of  Eq.  (1),  as  a 
function  of  the  relative  orientation  of  the  ground-  and 
excited-state  hyperfine  Hamiltonians.  The  arrow  at  the 
lower  left  points  to  the  coordinate  with  the  best  fit. 


somewhat  larger,  at  70  and  60  kHz,  respectively. 
The  results  for  the  'D2  excited-state  parameters  are 
summarized  in  Table  I. 

The  expression  for  echo  modulation  is  given 
analytically  by  Eq.  (1).  Once  Fourier  analysis  has 
given  the  splittings  and  linewidths,  the  only  remain¬ 
ing  parameters  define  the  relative  orientation  of  the 
ground-  and  excited-state  effective  hyperfine  Ham¬ 
iltonian  principal  axes.  This  Hamiltonian  may  be 
written 

^gui  =  ^g(e> 

(2) 


4(Zt  + 


Vgie> 


yin' 


The  upper  and  lower  case  subscripts  emphasize  that 
the  ground  (g)  and  excited  ie)  state  principal  axes 
need  not  be  the  same,  although  Neumann’s  princi¬ 
ple  requires  that  they  have  at  least  one  common  axis 
parallel  to  the  local  C2  symmetry  axis.  This  con¬ 
straint  significantly  reduces  the  number  of  possible 
orientations,  and  we  may  represent  these  possibili¬ 


TABLE  1.  Values  for  the  'D2  excited-state  hyperfine 
frequencies  and  nuclear  linewidths  obtained  from  the 
Fourier  transform  of  the  data  in  Fig.  1  The  linewidths 
have  been  corrected  for  the  effect  of  the  finite  truncation 
interval  (Ref.  7). 


Full  width  at 

Resonance 

Frequency 

half  maximum 

0)\ 

3.72+0.02  MHz 

70+20  kHz 

0)2 

4.79  +  0.02  MHz 

60+  20  kHz 

M\z 

8.51+0.02  MHz 

30+30  kHz 

ties  on  a  two-dimensional  map.7  The  coordinates  of 
such  a  map  give  the  particular  choice  of  common 
axes,  and  the  angle  0  through  which  one  set  of  axes 
is  rotated  with  respect  to  the  other  set  about  the 
common  axis.  Since  the  Hamiltonian  in  Eq.  (2)  is 
insensitive  to  the  sign  of  the  coordinates,  an  orienta¬ 
tion  of,  for  example,  x\\X  is  equivalent  to  x| |  —X. 
We  have  calculated  the  echo  modulation  curve,  Eq. 
(1),  in  increments  of  five  degrees  in  the  rotation  an¬ 
gle  0.  Figure  3  presents  a  map  where  for  each 
orientation  coordinate  a  box  is  drawn  whose  width 
is  proportional  to  the  minimum  mean-square  devia¬ 
tion  between  our  data  and  the  calculated  curves. 
Only  the  first  microsecond  of  data  was  used  in  the 
fitting  procedure  for  Fig.  3.  Beyond  1  /zsec,  the 
data  are  dominated  by  the  excited-state  hyperfine 
splittings  due  to  the  relative  linewidths  of  ground- 
and  excited-state  levels  and  including  these  data  in 
the  fit  considerably  reduces  the  selectivity  of  the 
map. 

Figure  4  shows  in  expanded  scale  the  first  three 
microseconds  of  (a)  our  data  and  (b)  the  theoretical 
curve  calculated  from  Eq.  (1)  using  the  orientation 
coordinate  giving  the  best  fit  to  the  data  (as  indicat¬ 
ed  by  an  arrow  in  the  lower  left  of  Fig.  3).  The 
curves  of  Fig.  4,  having  been  multiplied  by  the  fac¬ 
tor  exp(  +  4T/r2),  do  not  include  the  effect  of  the 
overall  optical  linewidth.  The  slight  rise  in  the  echo 
peaks  of  Fig.  4(a)  is  an  artifact  in  the  data  due  to 
detector  saturation  at  short  pulse  separation.  This 
saturation  may  also  be  seen  in  Fig.  1,  where  the  ini¬ 
tial  envelope  decay  is  somewhat  slower  than  the 
corresponding  decay  for  the  rest  of  the  data.  This 
artifact  introduces  some  uncertainty  in  the  estimate 
of  optical  linewidth  but  seems  to  have  negligible  ef¬ 
fect  on  the  Fourier  transform  and  overall  fit  since 
another  data  set  with  smaller  overall  signal-to-noise 
ratio,  but  without  the  saturation,  gives  similar  re¬ 
sults  to  those  presented  here. 

The  map  in  Fig.  3  shows  a  high  degree  of  sym¬ 
metry  and  concommitant  lack  of  selectivity  for 
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determining  the  orientation  parameter.  To  explore 
the  extent  to  which  this  symmetry  is  intrinsic  to  the 
form  of  Eq.  (1),  we  have  substituted  for  the  data 
used  to  calculate  the  map  of  Fig.  3  the  values  of  the 
theoretical  curve  corresponding  to  Fig.  4(b),  produc¬ 
ing  the  map  shown  in  Fig.  5.  The  scale  in  Fig.  5 
was  chosen  to  make  the  box  corresponding  to  the 
coordinate  with  the  worst  Fit  in  Fig.  5  approximate¬ 
ly  the  same  size  as  the  box  with  the  worst  fit  in  Fig. 
3.  It  is  apparent  that  the  high  degree  of  symmetry 
persists.  To  estimate  the  extent  to  which  our 
signal-to-noise  ratio  would  need  to  improve  and  so 
allow  us  to  select  one  orientation  with  certainty,  we 
replotted  Fig.  5  with  the  box  width  scale  expanded 
by  a  factor  of  10,  as  shown  in  Fig.  6.  It  appears 
that  an  improvement  in  the  signal-to-noise  ratio  by 
a  factor  of  10  or  20  would  have  made  possible  a 
more  definitive  statement  about  the  orientation. 
This  improvement  will  probably  come  from  the  use 
of  frequency-  and  amplitude-stabilized  dye  lasers. 
Frequency  jitter  in  particular  causes  considerable 
fluctuation  in  echo  intensity,  presently  averaged  by 
the  data-taking  system,  and  such  averaging  tends  to 
wash  out  some  of  the  fine  details  of  the  modulation 
pattern.  Despite  the  high  symmetry,  we  can  draw  a 
number  of  conclusions  about  the  relative  orienta- 
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tion.  All  configurations  with  all  axes  parallel  are 
ruled  out,  as  are  those  with  the  excited-state  Z  axis 
parallel  to  C2.  These  conclusions  are  similar  to  the 
measurement  by  Chen  et  al.1  of  the  relative  orienta¬ 
tion  of  the  3//4  and  3P0  states.  Since  the  g  tensor  is 
nearly  isotropic  in  the  'D2  state10  and  absent  com¬ 
pletely  in  the  *P0  state,  the  relative  orientation 
parameter  must  be  primarily  determined  by  the  re¬ 
lative  orientation  of  the  pure  quadrupole  tensor 
with  respect  to  the  ground-state  hyperfine  tensor. 

Our  photon  echo  modulation  analysis  of  the  'Z)2- 
3//4  Pr3+:LaF)  transition  hyperfine  structure  gives 
a  new  time-domain  measurement  of  the  splittings 
and  linewidths  which  accords  well  with  previous 
measurements.  We  also  present  new  information  on 
the  relative  orientation  of  the  effective  hyperfine 
Hamiltonians. 
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We  comment  on  the  spectral  modification  and  temporal  reshaping  of  a  narrow-band  light  pulse  propagating 
through  an  atomic  absorber  that  has  a  still  narrower  spectral  hole.  The  regime  of  propagation  narrowing  is 
distinguished  from  the  regime  of  Beer's  law. 
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I.  INTRODUCTION 

Hole  burning  is  the  term  for  the  selective  photo¬ 
excitation  of  atoms  or  molecules  in  a  small  por¬ 
tion  of  an  inhomogeneously  broadened  optical  ab¬ 
sorption  line.  In  this  region  of  the  line,  following 
the  photoexcitation,  there  will  be  fewer  than  the 
normal  number  of  atoms  or  molecules  in  their 
ground  states.  A  subsequent  narrow-band  weak 
probe  pulse  will  suffer  less  than  normal  attenua¬ 
tion  if  its  center  frequency  lies  within  the  excita¬ 
tion  region.  The  resulting  absorption  curve 
mapped  out  by  scanning  the  frequency  of  the  probe 
beam  has  a  hole  in  it.  Figure  1  shows  an  example. 
All  of  this  is  well  known. 1 

In  the  present  note  we  address  a  question  that 
comes  up  when  the  weak  probe  pulse  has  a  spec¬ 
tral  width  greater  than  the  width  of  the  hole,  as 
shown  in  Figure  2.  This  question  is:  What  is  the 
temporal  form  of  the  transmitted  pulse?  The  un¬ 
certainty  relation  A i/At  suggests  that  At  in¬ 


creases  because  At/  decreases,  but  in  order  to 
answer  the  question  properly  it  turns  out  to  be 
necessary  to  consider  the  propagation  of  the 
probe  pulse  through  a  finite  depth  of  absorber.  If 
the  probe  pulse  is  weak  enough  to  have  negligible 
effect  on  the  hole,  then  the  question  can  be  an¬ 
swered  analytically  in  closed  form. 

II.  PROPAGATION  OF  WEAK  OPTICAL  PULSES 
THROUGH  BROADBAND  INHOMOGENEOUSLY 
BROADENED  ABSORBERS 

In  the  familiar  two-level-atom  model  for  an 
absorber,2  transmission  and  absorption  of  light  are 
governed  by  the  equation 

(8/32  +  d/bct)S(t,z)  =  2  vi(w/c)<P(t,z) ,  (1) 

where  S  and  <P  are  the  complex  envelopes  of  the 
electric  field  strength  and  the  atomic  polarization 
density.  An  expression  for  (P  involves  the  off- 
diagonal  part  of  the  two-level  atom’s  density 


FIG.  1.  A  Doppler-broadened  absorption  line  with  a 
homogeneously  broadened  hole. 


FIG.  2.  The  initial  probe  pulse  spectrum  (shaded) 
overlaid  on  an  absorption  line  with  a  much  narrower 
spectral  hole. 
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matrix  in  the  rotating  wave  approximation,  averaged 
over  the  Doppler  distribution  of  atomic  veloci¬ 
ties: 

<P  (t,z)  =  Nd(rl2(t,z;  i’))„ .  (2) 

Here,  N  is  the  density  of  atoms  and  d  is  the  mag¬ 
nitude  of  the  atomic  transition  dipole  matrix  ele¬ 
ment  along  the  direction  of  the  electric  field  vec¬ 
tor. 

If,  as  we  assume,  the  electric  field  of  the  probe 
pulse  is  too  weak  to  alter  the  shape  of  the  inhomo¬ 
geneous  line  (with  the  hole  already  burned  into  it), 
then  the  equation  obeyed  by  rl2it,z;  v)  is  simply 

drl2/at  =  -  (0  +  *a)»T2  +  (i/K)dS.  (3) 

Here,  A  =  A(u)  is  the  detuning  between  pulse  fre¬ 
quency  and  the  Doppler -shifted  atomic  resonance 
frequency  and  0  is  the  halfwidth  of  the  homogeneous 
component  of  the  atomic  absorption  line.  An 
average  over  A ’s  is  equivalent  to  an  average  over 
v’s. 

Because  Eqs.  (1)  and  (3)  are  linear  they  are 
trivially  solved  by  Fourier  transform  methods. 

We  define 

Sit, 2)=  f  e(v,z)e~ivtdv/2n ,  (4a) 

rl2(t,z)=  f  p(v,z)e~ivtdv/2v ,  (4b) 

and  find  the  equation 

[d/dz-iv/c  -if(v)+  \  g(v)]e(v,z)  =  0,  (5) 


linewidth  is  appreciably  greater  than  6v,„  we  can 
take  g{v)  =  g{ 0)  for  all  significant  frequencies 
( I  v  |  •-  5v0).  Then  the  only  significant  frequency 
dependence  is  due  to  the  Gaussian  factor  that  came 
from  the  incident  pulse.  Thus  the  pulse  retains 
its  Gaussian  shape  and  is  simply  attenuated  uni¬ 
formly  at  all  frequencies  according  Beer’s  law: 

|  e(v,z) \2=  |  S0\2e~,v /6t,° > 2 e~‘(0 ** .  (10) 

All  this  is  well  known,  and  Crisp3  has  given 
interesting  examples  showing  dramatic  departures 
from  the  kind  of  behavior  implied  by  (10)  when 
the  pulse  bandwidth  is  much  greater,  rather  than 
less,  than  the  Doppler  linewidth. 

Hi.  PROPAGATION  THROUGH  A  SPECTRAL  HOLE 

With  the  background  formulas  derived  in  Sec.  II, 
it  is  easy  to  take  into  account  a  narrow  hole  in  the 
Doppler  line.  For  simplicity  we  will  take  the 
shape  of  the  Doppler  curve,  as  well  as  the  hole, 
to  be  Lorentzian.  This  is  realistic  as  far  as  the 
hole  is  concerned  and  not  a  bad  approximation 
for  the  Doppler  curve  because  we  will  be  con¬ 
cerned  only  with  its  center  portion.  Another 
simplication  will  be  to  take  the  centers  of  the 
spectral  curves  of  the  incident  pulse,  the  Dop¬ 
pler  line,  and  the  hole  all  to  coincide  at  v  =  0,  as 
shown  in  Fig.  2.  Under  these  conditions  the 
angular  bracket  in  Eqs.  (6)  is  to  be  interpreted  as 
an  average  over  Doppler  detuning  A  with  the 
normalized  weight  function 


where  /  and  g  are  the  dispersive  and  absorptive 
parts  of  the  dipole  reaction  field: 


/M-  2GIm(^+  ;(*_„)), 

(6a) 

2gM~  2GRe(|3+i(A-K))  ’ 

(6b) 

Here,  G  is  the  primitive  attenuation  parameter  for 
the  problem 

G=  4it Nd2u/Hc  ,  (7) 

and  the  angular  brackets  denote  the  Doppler 
average. 

The  solution  to  Eq.  (5)  is  immediate: 

e{v,z)=  S0e-l,2«',6v  0>2  w  (8) 

We  have  assumed  that  the  incident  pulse  is  Gaus¬ 
sian  in  time,  with  bandwidth  Si/0  and  temporal 
length  2  v/6i>0. 

The  pulse  Fourier  energy  spectrum  obviously 
changes  with  propagation: 

\e(v,z)\l=  \  S0V  ev,,,v<>'2  e"”* .  (9) 

However,  in  the  usual  case,  when  the  Doppler 


pi A)  = 


0*  1 
u  Az+0+)2’ 


where  0*  is  obviously  now  the  Doppler  width.  Then 
we  can  evaluate  the  brackets  in  Eqs.  (6)  by  explicit 
integration  and  find 


f(v)= 

5g(F)=  \a  (- 


v/0*  _  v/Ph 

1  +  (u/0*)2  1  +  (v/0H) 


1  _  1 

1  +  iv/P*)2  1  +  (v/ft,)2 


)■ 


(Ha) 

(lib) 


where  a  is  the  normal  Doppler  absorption  coeffi¬ 
cient  at  line  center: 


a  =  4  w  Nd2w  /ftc  0*  ,  (12) 

and  0*  and  0B  are  the  half -widths  of  the  Dopper 
line  and  the  hole,  respectively.  We  are  principally 
interested  in  the  limits 


0*»  6v0»  0y»  0 .  (13) 

In  other  words,  we  assume  the  Doppler  line  to 
be  much  broader  than  any  other  spectral  feature, 
and  the  width  of  the  incident  pulse  to  be  much 
broader  than  the  hole  in  the  Doppler  line.  The 
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underlying  homogeneous  width  is  taken  to  be 
negligibly  small  compared  to  all  the  other  widths. 
These  conditions  are  qualitatively  met  by  the 
curves  in  Fig.  2. 

Given  relation  (13),  there  are  two  distinct  re¬ 
gions  of  the  spectrum: 

(a)  The  region  of  propagation  narrowing,  where 

v  ^  Pg. 

(b)  The  region  erf  Beer's  law,  where  v»  &H. 

We  have  discussed  region  (b)  briefly  in  Sec.  II. 
Region  (a)  is  characterized  by  a  simple  approxi¬ 
mate  expression  for  g(u): 

g(v)=*a(  v/$gf.  (14) 

Thus,  in  region  (a)  we  have 

|e(v,*)|2=  (15) 

where 

6vU)=  pB/[az  +  (|3tf/6i/0)2]1/2 .  (16) 

That  is,  in  the  region  of  propagation  narrowing 
the  spectrum  is  Gaussian  but  with  a  variable 
width.  Note  that  at  atz  =  0  the  width  is  equal  to 
the  initial  pulse  width  6i/0  and  only  becomes  small¬ 
er  due  to  propagation.  The  decrease  can,  however, 
be  rapid  because  PH/5v0  can  be  much  less  than 
unity.  The  threshold  for  propagation  narrowing 
occurs  at  az~  <&„/ 5v0)2. 

In  the  region  of  Beer’s  law  the  spectrum  is 
also  Gaussian,  but  with  width  6 u0,  independent 
of  atz. 

To  interpolate  between  regions  (a)  and  (b)  re¬ 
quires  the  full  expression  (lib)  for  g(v).  In  Fig. 

3  we  show  the  full  \e(v,z)\2  as  a  function  of  i>  for 
several  values  of  atz,  given  the  ratios  p*/6v0 
=  /TO  and  pa/6 v0=  1/10.  The  two  distinct  Gaus- 
sians  appropriate  to  spectral  regions  (a)  and  (b) 
are  evident.  Figure  4  shows  the  effect  of  the  hole 


v/tv0 

FIG.  4.  The  spectrum  of  the  probe  pulse,  as  tu  Fig. 

3 ,  except  that  the  normal  component  of  absorption  has 
been  removed. 

alone  by  plotting  \e(v,z)\2  without  the  first  term 
of  g(v).  In  effect,  Fig.  4  shows  that  the  relative 
importance  of  the  hole  grows  with  propagation, 
but,  as  expected,  only  near  line  center.  The  broad 
base  of  all  erf  the  curves  is  simply  exp[-  (v/6i/0)2]. 

IV.  TIME-DEPENDENT  FEATURES 

The  exact  expressions  obtained  in  Sec.  Ill  for 
\e(v,z)\2,  as  well  as  the  approximate  relations 
valid  in  the  regions  of  propagation  narrowing  and 
Beer’s  law,  permit  good  qualitative  estimates  of 
the  time  dependence  of  the  transmitted  pulse. 
However,  |e(v,z)|2  does  not  provide  everything, 
because  the  consequences  oif(u)  are  not  included. 
In  this  section  we  present  the  results  of  numerical 
computation  of  the  Fourier  transform  (4a),  there¬ 
by  giving  the  full  space -time  behavior  of  S(t,z). 

Figure  5  shows  the  temporal  behavior  of  the 
transmitted  pulse  after  propagation  to  az  =  2.  A 
clear  departure  from  the  Gaussian  shape  of  the 


FIG.  3.  The  spectrum  of  the  probe  pulse  [eiv.z)^  at 
a  succession  of  propagation  depths  az«o.O,  0.5,  l.o, 
and  2.0.  The  horizontal  axis  is  in  units  of  5v#. 


FIG.  5.  The  probe  pulse  intensity  |5«,z)|2  in  arbi¬ 
trary  units  as  a  function  of  time,  after  propagating  two 
absorption  depths  into  the  medium.  The  horizontal  axis 
is  in  units  of  (Ai0)'* . 
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FIG.  6.  The  probe  pulse  Intensity  )<S(/ , *) I1  as  a  func¬ 
tion  of  time  at  two  different  spatial  positions,  at  =0  and 
2.  The  vertical  axis  is  logarithmic,  but  to  the  same 
scale  as  in  Fig.  5,  and  the  units  of  time  are  (6 v0) "* . 

input  pulse  is  evident  in  the  trailing  portions  of 
the  transmitted  pulse.  In  Fig.  6  we  show  the  data 
of  Fig.  5  again,  but  with  the  z  =  0  curve  added  for 
comparison.  At  bvj  -  2  the  trailing  edge  of  the 
transmitted  pulse  becomes  stronger  than  the 
tail  of  the  incident  pulse,  although  it  is  two  orders 
of  magnitude  weaker  than  the  incident  pulse  peak. 

Finally,  Fig.  7  shows  the  transmitted  pulse  at 
az  =  6.  It  is  seen  that  the  bulk  of  the  pulse  energy 
is  now  in  the  long  tail,  and  one  can  begin  to  say 
that  a  single  frequency -time  uncertainty  relation 
again  describes  the  pulse  reasonably  well.  The 
temporal  width  of  the  tail  is  roughly  V"3  times  the 
width  shown  in  Fig.  5,  in  accord  with  (16)  since 
(Ph/6v0)2  is  negligible.  Moreover,  we  note  that 
the  resonant  interaction  of  pulse  and  atoms,  during 
the  pulse  transmission,  is  fully  coherent  in  our 
model.  This  is  responsible  for  the  slight  inter¬ 
ference  minimum  that  occurs  between  the  two 
components  of  the  pulse. 

V.  DISCUSSION 

The  spectral  and  temporal  changes  predicted 
by  solution  (8)  are  illustrated  in  the  figures. 

They  are  a  consequence  of  purely  linear  disper¬ 
sion  theory  and  are  novel  only  in  the  sense  that 
classic  discussions4  of  light  propagation  in  dis¬ 
persive  media  do  not  appear  to  have  included 
treatments  either  of  purely  inhomogeneous  line 
broadening  or  of  lines  with  holes  in  them.  Ex¬ 
perimental  observations  of  the  predicted  peak 
delay  and  pulse  lenghtening  do  not  appear  to  have 
been  reported  either. 

Another  view  of  these  results  is  obtained  by  re¬ 
garding  our  model  of  an  absorption  line  with  a 
hole  as  a  continuous -band  interference  filter.  In 
one  sense  it  is  a  pure  interference  filter;  no 


FIG.  7.  The  probe  pulse  intensity  |5(/,*)|*  as  a  func¬ 
tion  of  time  after  propagating  six  absorption  depths  into 
the  medium.  The  vertical  axis  is  to  the  same  scale  as 
Fig.  5,  and  the  units  of  time  are  (fir,)"'. 

energy  loss  mechanisms  are  included  in  the  model. 
However,  conversion  of  pulse  energy  into  loss- 
free  atomic  dipole  oscillations  all  across  the  Dop¬ 
pler  line  will  nevertheless  lead  to  pulse  decay, 
and  it  is  more  realistic  to  speak  of  an  inter¬ 
ference  filter,  as  opposed  to  an  absorption  filter, 
when  the  effect  of  f(v)  is  much  more  important 
than  that  of  g(v).  Formulas  (11)  show  that  for 
v  «  this  condition  is  well  satisfied: 

/O')*  I  an'//3„,  (17a) 

IgO*)*!  (17b) 

This  is  a  consequence  of  the  hole  in  the  line,  of 
course.  For  a  normal  line  without  a  hole,  one 
has  g»  f  for  all  frequencies  near  the  center. 

Finally,  because  of  the  linearity  of  the  model, 
it  is  also  possible  to  treat  the  pulse  as  if  it  were 
interacting  with  two  entirely  separate  groups  of 
atoms.  The  first  group  comprises  the  usual  di¬ 
pole  oscillators  in  the  full  Doppler  line,  and  the 
second  group  is  a  set  of  “negative”  oscillators 
occupying  a  region  of  width  0H  at  line  center.  The 
negative  oscillators  emit  rather  than  absorb  light. 
It  is  the  emission  of  the  "negative”  oscillators 
that  causes  the  growth,  with  increasing  propaga¬ 
tion  distance,  of  the  peaks  of  the  curves  in  Fig.  4. 
In  this  view  of  the  model  it  is  the  interference 
between  the  positive  and  negative  oscillators  at 
line  center  that  causes  the  pulse  lengthening.  In 
classical  theory  the  negative  dipoles  can  never 
be  more  “negative”  than  the  real  dipoles  are 
“positive,”  because  they  are  designed  only  to 
cancel  the  real  dipoles  in  a  certain  spectral  re¬ 
gion.  In  the  quantum  theory  the  negative  atoms 
are  not  so  severely  restricted.  For  example, 
the  hole  has  been  constructed  in  our  treatment 
so  that  g( 0)  =  0.  That  is,  we  have  assumed  the 
absence  of  absorbing  dipoles  at  line  center.  How- 
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ever,  the  original  preparation  of  the  hole  could 
have  been  arranged  to  include  a  degree  of  atomic 
inversion  at  line  center.  In  this  case  g(0)  <  0,  and 
gain  rather  than  loss  would  be  expected.  The 
“negative”  oscillators  in  the  model  would  then 
more  than  cancel  the  positive  ones.  We  will  not 
explore  this  case. 

ACKNOWLEDGMENTS 

This  research  was  partially  supported  by  the 
Department  of  Energy,  Contract  No.  DE-AC03- 


'See,  for  example,  discussions  in  A.  Yariv,  Quantum 
Electronics  (Wiley,  New  York,  1967),  Chap.  8; 

M.  Sargent  HI,  M.  O.  Scully,  and  W.  E.  Lamb,  Jr., 
Laser  Physics  (Addison-Wesley,  Reading,  Mass., 
1974),  Chap.  10;  and  L.  Allen  and  J.  H.  Eberly,  Opti¬ 
cal  Resonance  and  Two-Level  Atoms  (Wiley,  New 
York,  1975),  Chap.  6. 


76DP01118,  the  Joint  Services  Electronics  Pro¬ 
gram,  Contract  No.  DAAG29-79-C-0079,  the 
National  Research  Council  of  Canada,  the  National 
Science  Foundation  and  Army  Research  Office, 
Grant  No.  NSF-DMR80-06966,  and  the  Office  of 
Naval  Research,  Contract  No.  N00014-79-0666. 

We  are  grateful  to  J.  J.  Yeh  and  B.  J.  Herman  for 
assistance  with  the  numerical  integration  required 
for  Figs.  5-7,  and  we  acknowledge  conversations 
with  M.  A.  Johnson,  B.  W.  Shore,  S.  Stenholm, 

A.  Szdke,  and  J.  J.  Yeh. 


2We  follow  the  conventions  of  L.  Allen  and  J.  H.  Eberly, 
see  Ref.  1. 

3M.  D.  Crisp,  Phys.  Rev.  A  1,  1604  (1970).  See  also 
Ref.  2,  Sec.  1.8. 

4A.  Sommerfeld,  Optics  (Academic,  New  York,  1972); 
L.  Brillouin,  Wave  Propagation  and  Group  Velocity 
(Academic,  New  York,  1960). 


\-  — .  I-—.  4. 


